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I. INTRODUCTION 



The goal of this work is to find the effective permittivity tensor for an arbitrary doubly 
periodic array of circular cylinders when dielectric properties of the cylinders and the host 
medium are complex-valued. The problem of determining effective properties of a periodic 
array of inclusions was first considered by Rayleigh 1 . His method was extended 2 ' 3 for reg- 
ular arrays of circular cylinders. A method of functional equations 4,5 employing analytic 
functions was used to find an expression of the permittivity tensor for small volume fraction 
of inclusions. For rectangular lattice of inclusions an efficient method based on the use of 
elliptic functions was suggested in Ref. 6. Numerical evaluation of complex- valued effective 
permittivity can be found in 7 , 8 as well as in 9 where one can find the bounds on the transport 
coefficients of a two-dimensional composite. 

The solution of the problem in question consists of two steps. First, we construct a 
quasiperiodic potential using a combination of the Weierstrass ^-function and its derivatives. 
This ensures periodicity of the electric field in the whole plane and, as a result, avoids the 
problem of the summation of a conditionally convergent series. This approach is similar 
to that in Ref. 10 applied to biharmonic problems of the theory of elasticity. Second, we 
determine the average electric field and electric displacement within one parallelogram of 
periods and find an exact expression of the effective permittivity tensor that relates the two 
quantities. This paper is an extension of the author's previous article 11 to the case when 
physical properties of the component constituents are complex-valued. The use of complex 
quantities allows us to greatly simplify the solution of the problem. 

II. DERIVATION OF PERIODIC POTENTIAL 

Suppose that a periodic lattice of identical circular inclusions of radius a is introduced 
into a uniform complex electric field E applied in the plane perpendicular to the cylinder 
axes. The nodes of the lattice in the complex plane are generated by a pair of vectors 2ri 
and 2r 2 , Im ^ > (see Figure 1). In polar coordinates the complex- valued potential u{r, •&) 
has the following properties: 




u in in the inclusion, 



u ex in the medium, 



(1) 
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where e(r,u) = e'{r,uj) + ie"{r,uS) is the tensor of complex permittivity at frequency 
u. Hereafter prime and double prime denote the real and imaginary part of the quantity, 
respectively. 

Potentials Ui n and u ex cannot be both analytic, and therefore we represent u in the form 

■ z \ 2n+l 



u in {z) = Ea\^ A n[-J +B n i K 

n=0 L 



~ x 2n+l 



Uex(z) = -Ez + Ea^2 



x> a 2n+1 



n=0 



(2n)! 



\z\ < a, 
[C n ( (2n \z) + D n ^ 2n \z)] . 



\z\ > a, 



(2) 
(3) 



where A n , B n , C n , and D n are unknown complex dimensionless coefficients, z stands for 
the complex conjugation, and ((z) is Weierstrass' ^-function 12 



m,n 



1 



1 : 
+ 



z - P P ' P 2 

" 1 m,n 1 m,n 1 m n 



(4) 



Here C,^ 2n \z) denotes derivative of order 2n, and P m ,n = 2mri + 2riT2. Prime in the sum 
means that summation is extended over all pairs m, n except m = n = 0. 
On the boundary r = a of inclusion we impose continuity conditions 



u] = 0, 

du 
dr 

where brackets [•] denote the jump of the enclosed quantity across the interface. 



(5) 
(6) 




FIG. 1. Circular inclusions of radius a arranged in a periodic lattice with periods 2t\ and 2t2- 



Below we will use some properties of ^-function 12 ' 13 . It has the quasiperiodicity property 



C(z + 2r k ) = C(z) + 2 Vk , Vk = ((T k ), k — 1,2, 



(7) 



where constants r]i and r] 2 are related by the Legendre identity 

71% 

Vin -V2T1 = —. (8) 

The derivative of ((z) is a periodic function and is expressed through Weierstrass elliptic 
function p{z) by 

C(z) = -p(z). (9) 

This property ensures the electric field to be periodic in the medium, while (7) guarantee 
that the potential changes by a constant value in the direction of either 2t\ or 2r 2 . 

To satisfy conditions (5)- (6) on the inclusion surface we expand ((z) and its even deriva- 
tives in a Laurent series 

C^(z)- (2n)l V* (2n + 2fc + l)! 

fc=0 

where 

s * = £pt> fe = 2,3,.... (11) 

Sums (11) contain only even powers of P m , n since for every point P m , n = 2mri + 2nr 2 on 
the lattice there exists symmetric point —P m ,n and the sums with odd powers vanish. Also, 
if the periods 2ri and 2r 2 of the lattice are fixed, sums Sk remain bounded as k — > 00. 
We introduce new variables 

On = A n + i? n , b n = A n — B n , (12) 
c n = C n + D n , d n = C n — D n . (13) 

In the new variables, compliance with conditions (5)- (6) leads to a system of equations 

If If 

E(2n + 2m+l)! 2 „ +2m+2 r / 1r > 

77^ , T\TJr> M « Sn+m+l Cm = «0 n , , (15) 

m=0 (2n + l)!(2m)! 

d n + ajr ^ + ^1 a 2n+2m+2 s„+ m+ i d m = -a5 nfi . (16) 

^ (2n + l)!(2m)! 



where 



« = — — — • (17) 
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Using notation 



1 = 


min 




a 


h = 


T 


Sk = 


£' 




n.m 



1 



p 

both equations (15)-(16) then can be written as 



r2|,2|7i-T 2 |}, (18) 

(19) 

, fc = 2,3,..., S 1 = 0, (20) 



x„ - J] x m /i 2 " +2m+2 = ay 5 n , , (21) 

m=0 

where 

_ (2n + 2m + l)! _ _ 

(zn + l)!(2mj! 

The system (21) has the same structure as that considered in 11 . Therefore, its unique solution 
exists for sufficiently small values of h and can be obtained in the form of a convergent power 
series 

oo 

x n = ay 5 n , + Pn, m h 2n+2m+2 , (23) 

m=0 

where 

p nfi = aG nfi y, 

Pn,k = J2 Gn 

k = l,2,.... (24) 

m=0 

Calculation of Co and do from (23)-(24) gives 



c = aX, d Q = —a/i, (25) 



where 



,16 



A = 1 + 3a 2 S 2 2 h 8 + 5a 2 Sjh 12 + 30a 3 5 2 5 3 /i 14 + (9a 4 S 2 4 + 7a 2 Sj) h 1 
+ 210a 3 S 2 S 3 S 4 h 18 + (330a 4 S 2 S 2 + 9a 2 S 2 )h 20 + 0(h 22 ), (26) 

H = 1 + 3a 2 S 2 2 /i 8 + 5a 2 S 3 V 2 - 30a 3 5 2 2 S 3 /i 14 + (9a 4 S 2 4 + 7a 2 Sl) h 16 
- 210a 3 S 2 S 3 S 4 h 18 + (330a 4 S 2 S 2 + 9a 2 S 2 )h 20 + 0(h 22 ). (27) 
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III. CALCULATION OF THE EFFECTIVE PERMITTIVITY TENSOR 



Effective permittivity tensor e* relates the average electric displacement (£>) and the 
average electric field (E) 

(D)=e*(E). (28) 

Observe that 



(D} = - DdS=^ . 

D Js D J Si. 



E m dS + / E ex dS. 

** Sex 



Calculation of the above integrals gives 11 

Ej, n dS 



2na 2 e f 



S J S 

1 



ME, 



E er dS = E — &ME 



where 





A 







1 


M = 


, E = E 







n 




% 



and matrix *P depends only on geometry of the lattice 

Re77ilm2r 2 — Im?7 1 Im2T2 
— Im r/ilm 2r 2 7r — Re r^Im 2r 2 

Computing the average electric field (E) we obtain 



(E) 



I - 



2aa 2 



&M] E. 



S J 

From here and (30)-(31) we derive the effective permittivity tensor 

e* = e ex {I + n5M (I - 5SPMy l } , 

where 

2aa 2 
o = — - — . 



(29) 

(30) 
(31) 

(32) 



(33) 



(34) 



(35) 



(36) 



Explicit computation of the 2x2 inverse matrix in (35) shows that e* is symmetric. If 
\\5\PM || < 1 then e* can be expanded in a convergent series 



e* = e ex I + e ex n6M (5&M) n . 



(37) 



n=0 



In what follows we use expression (35) for calculation of the effective permittivity tensor for 
specific lattices. 
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IV. EXAMPLES 



A. Square lattice 

For the square lattice we put 2ti = I, 2r 2 = it (see Figure 2). Then one can find 13 that 

7T 7T 

?7i = — . From (33) we obtain & = —I. All lattice sums (20) are real with the only non-zero 













2x 2 











FIG. 2. Square lattice of inclusions of radii a with periods 2t\ = I and 2r 2 = it. 

being S^fcj A; = 1,2,.... Calculation of matrix M in (32) gives X — /i. Hence tensor e* = e*I 
is isotropic with 

1 + aA/ 



e = e f 



1-aXf 



(38) 



where 



A = 1 + 3a 2 Sih 8 + a 2 (9a 2 S% + 7S 2 A ) h 16 + (27a 6 S% + 21a 4 S 2 2 S 4 2 + 2205a 4 S 2 2 S 4 2 
+ 21a 4 S 2 2 Sl + lla 2 S 2 6 ) h 2A + 0(h 32 ). (39) 

Here f = irh 2 is the volume fraction of the inclusions, h — —, So — > -, rr = 3.15121, 

n,m 

S 4 = V'- — — = 4.25577, S 6 = V' 7 ^— -r^ = 3.93885 correct to five decimal 

(m + mr (m + zn) u 

n,m n,m 

places. Dependence of the real and imaginary parts of e* on the volume fraction of the 
cylinders is shown in Figures 3 and 4 for e in = 30 — 0.3i, e ex = 1 — hi and e in = 1 — 8i, 
e ex = 2 — 0.3i, respectively. In both examples the dielectric constants are chosen in such 
a way to make either real or imaginary part of e* nonmonotonic. Formula (38) gives an 



excellent agreement between numerical evaluation of e* form (15), (25) and calculation by 
(39). 




(a) (b) 

FIG. 3. Dependence of the real (a) and imaginary (b) part of the complex effective dielectric 
constant e* on the volume fraction / evaluated numerically (solid blue line) and by formula (38) 
(red circles) for a square array of cylinders. We used. Ei n = 30 — 0.3i and e ex = 1 — hi. 
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(b) 



FIG. 4. Dependence of the real (a) and imaginary (b) part of the complex effective dielectric 
constant e* on the volume fraction / evaluated numerically (solid blue line) and by formula (38) 
(red circles) for a square array of cylinders when £j n = 1 — 8? and e ex = 2 — 0.3i. 
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B. Regular triangular lattice 



For triangular lattice (see Figure 5) we put 2ti = I, 2r 2 
constant r\\ 

e* = e*I. We can use (38) with 



— ' K —=. As a result, matrices \P = —I and M 



£e ni/3 . Then we find 13 the 
XI are isotropic as well as 



A = 1 + 5a 2 S 2 h L2 + a 2 (25a 2 Sl + HS 2 6 )h 24 + 0(^ e ), 



(40) 



where / = -^=h 2 , h — ^, while the only non-zero real lattice sums are S 3 k, k = 1, 2, 

Here S 3 = V' — = 5.86303, S 6 = V' — — jg = 6.00964 correct to five 



n,m I m + ne 3 I 



n,m + ne 3* j 




FIG. 5. Regular triangular lattice of inclusions of radii a with periods 2t\ = £ and 2r 2 = ie 11 *^ . 

decimal places. Dependence of the real and imaginary parts of e* on the volume fraction of 
the cylinders is shown in Figures 6 and 7 for e in = 30 — 0.3i, e ex = 1 — hi and Ei n = 1 — 8i, 
e ex = 2 — 0.3i, respectively. In both examples the dielectric constants are chosen in such 
a way to make either real or imaginary part of e* nonmonotonic. Formula (38) gives an 
excellent agreement between numerical evaluation of e* form (15), (25) and calculation by 
(39). 

C. Anisotropic lattice 



As an example of anisotropic lattice we consider the latice with the periods 2r x = £, 
2r 2 = le 3m / 8 (see Figure 9). Then the constant % can be found numerically r/i = C( T i) ~ 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 14 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0^9 1 

f f 

(a) (b) 

FIG. 6. Dependence of the real (a) and imaginary (b) part of the complex effective dielectric con- 
stant e* of a regular triangular array of cylinders on their volume fraction / evaluated numerically 
(solid blue line) and by formula (38) (red circles) for e^ n = 30 — 0.3i and e ex = 1 — 5i. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

f f 



(a) (b) 

FIG. 7. Dependence of the real (a) and imaginary (b) part of the complex effective dielectric con- 
stant e* of a regular triangular array of cylinders on their volume fraction / evaluated numerically 
(solid blue line) and by formula (38) (red circles) for £i n = 1 — 8i and e ex = 2 — 0.3i. 



(1.73289 - 0.07888i)/£, and matrix <F in (33) becomes 



3.20196 0.14575 
0.14575 -0.06037 



(41) 
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Lattice sums for approximation of A and \i in (26)-(27) have the values 5(2) = 




1 



m + 
1 



ne s 



r 



1.01011+1. OlOlli, 5(3) = 




8 



0.87457i, 5(5) = 



n,m 



m + 
1 



n,m (m + ne 3 ™ 



) 
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of the effective dielectric tensor from (35) using the solution of (15)-(16) for numerical 
evaluation of matrix M and its approximation by (26)- (27). Figure (9) shows dependence 
of the components of the real and imaginary parts of e* on the volume fraction / of the 
cylinders. Lines denote component of the effective tensor obtained from numerical solution 
of (15)— (16), while symbols correspond to their approximation by (26)-(27). As in the 
previous cases, agreement between numerical evaluation and analytic approximation of the 
effective dielectric tensor is excellent even for high values of /. 



FIG. 8. Anisotropic lattice of inclusions of radii a with periods 2t\ = i and 2t2 = £e 3nt ^ 8 . 



We determine the effective complex permittivity of a two-dimensional composite, con- 
sisting of an arbitrary doubly periodic array of identical circular cylinders in a homogeneous 
matrix, and whose dielectric properties are complex-valued. Efficient formulas are provided 
to determine the effective complex permittivity tensor which are in excellent agreement 
with numerical calculations in the whole range of cylinder volume fraction variation. We 





V. CONCLUSION 



11 




(a) (b) 

FIG. 9. Dependence of the real (a) and imaginary (b) part of the complex effective dielectric tensor 
e* of a lattice with periods 2t\ = £, 2r 2 = £e 37 ™/ 8 on the volume fraction / of the cylinders. The 
solid blue line is e* l5 the dashed black line corresponds to £22 > the dotted red line is e* 2 evaluated 
numerically. Blue circles, black triangles, and red squares correspond to their approximation by 
(26)-(27), respectively. The dielectric constants of the materials are £i n = 30 — 0.3i and e ex = 1— hi. 

also show that in contrast to the real-valued case, the real and imaginary parts of the ef- 
fective complex-valued tensor can exhibit non-monotonic behavior as functions of volume 
fraction of cylinders, and can be either greater or less than that of the constituents. This 
raises the question about optimal design of such composites. 
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